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NON- STATION ARY INFINITE SERVER MODELS AND THEIR RELATIVES 



by 



Donald P. Gaver 

Naval Postgraduate School and 
Monterey, CA 93940 



John P. Lehoczky 
Carnegie-Mel Ion University 
Pittsburgh, PA 15213 



1 . Introduction 

The waiting-time model characterized by (i) Poisson arrivals 
to (ii) an unlimited number of servers, these characterized by 
(iii) independent service times of arbitrary distribution — usually 
called the M/G/°° system — has a special significance and utility 
that stems from the simplicity of its solution. That is, if N(t), 
t 0 , denotes the number of arrivals being served at time t, con- 
veniently referred to as system occupancy , and if N(0) = 0, then 
N(t) itself has the Poisson distribution. This fact is well-known 
when the Poisson arrival rate is a constant. A, and F(x) is the 
distribution of service times perhaps with finite mean E[S] = y ^ < co ; 
in this latter case the limiting distribution (t -*■ °°) is always 
Poisson with parameter E[S] = A/y; see Parzen (1962). Such a 
model approximately characterizes 

a) the number of occupied channels in a system of parallel, 
lightly-loaded telephone or communication channels (see 
Feller (1967) ) . 

b) the number of items undergoing repair in a logistics system, 

c) the number of vehicles simultaneously using a city's free- 
way system (see Newell (1966)). 

d) the number of drug "particles" inhabiting a particular organ 
(compartment) in a human or animal body at time t (see 
Gaver and Lehoczky (1977)). 



Research partially supported by National Science Foundation Grant 
Number MCS-77-07587 at the Naval Postgraduate School. 
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In this paper it is shown that N(t) possesses a 
Poisson distribution even if the arrival rate is time dependent, 
denoted by A(t) , and if the distribution of the service time 
S , of an arrival at t is itself dependent upon t, i.e. 
is F(x;t), t being a parameter. Furthermore, if arrivals 
occur in bunches — is compound Poisson — then N(t) itself 
has the compound Poisson distribution. Incidentally, if 
A ( t ) = A, a constant, arriving bunches are geometrically 
distributed, and service times are exponential, the limiting 
distribution of N(t) is shown to be the negative binomial. 

The simple methods used are extended to study multivariate 
process models, and also to shot-noise, storage, and Zipf's 
law. 



2 . A Backward Equation for System Occupancy: Poisson Arrivals 

The distribution of N (t) may be approached as follows. 

Let L\I(t,u) denote the system occupancy at t that arrived after 
time u, 0 < u < t; clearly one takes N(t,0) = N(t) , and 
N ( t , t ) = 0 with probability one. Let 

P n (t,u) = P{N(t,u) = n} (2.1) 

Hold t fixed, and consider the possible events that may occur in 
a time period of length du. Write F(x,t) = P{S t > x} = 1 - F(x,t) , 
and observe that either (a) no arrival occurs that remains a 
system occupant by time t, an event of probability 
1 - A(u)du F(t-u,u) + o(du), or (b) one arrival occurs and 
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remains an occupant at t, an event of probability 
A(u)du F(t-u,u) + o(du). All other possibilities are negligible. 
Thus 



p n (t,u + du) 



[1 - A(u)du F(t-u,u)] 
+ A (u) du F ( t-u , u) 



P n (t^u) 

P n _i (t , u) + o (du) 



( 2 . 2 ) 



Now subtract p (t,u) from each side, divide by du , and let 



du 0 . The result is the backward equation 



-dp n (t ,u) 
du 



-A (u) F (t-u, u) p n (t,u) + A(u) F (t-u , u) p n _^(t,u). 



n = 1,2,3, 



(2.3) 



and 



-dp Q (t,u) 
du 



■A (u) F (t-u , u) Pq ( t , u ) 



(2.4) 



The equation (2.4) may be integrated from u to t to give 



Pq ( t , u ) = exp [- / A ( v) F (t-v) dv] , 

u 



(2.5) 



and when u 0 this shows that 



p n (t,0) = P{N(t)=0} = exp[- / A (v) F(t-v,v)dv] . 

0 
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Invocation of the earlier equation (2.3) shows inductively that 



p n (t,0) = P{N(t) = n} 

t t _ 

= exp[- / A(v) F(t-v,v)dv] — r [-/ A(v) F(t-v,v)dv] n 
0 0 



( 2 . 6 ) 



and thus N(t) has a Poisson distribution. 

An approach via generating functions is easy and direct. 

Define 



. . . oo 

E[z (t,U) ] = g ( z , t ; u) = £ z n P{N(t,u) = n} . (2.7) 

n=0 

Note that the contribution to occupancy at time t from arrivals 
in (u, u + du) , AN(t,u) = N(t,u) - N(t, u + du) is independent 
of the contribution from arrivals after u + du by Poisson 
process properties, so the convolution property of the generating 
function leads to writing 

g ( z , t ; u ) 

= (z A(u)du F(t-u,u) + 1 - A(u)du F(t-u,u)] g(z,t; u + du) + o(du 

( 2 . 8 ) 



Next subtract g(z,t,u + du) from each side, divide by du 
and let du 0 . The result is 



dg(z,t;u) 

du 



(z-1) A(u)du F(t-u, u) g(z,t;u) 



(2.9) 



4 



This equation is immediately solved to produce 



t 

g(z,t;u) = exp[(z-l) / A (v) F(t-v,v) dv] , (2.10) 

u 

which is recognized to be the generating function of the Poisson 
distribution (2.6). 



3 . Compound Poisson (Bunched) Arrivals 

Next consider the situation in which the number of arrivals 
that occur together is random. In other words, bunches of 
arrivals occur together, and the bunch sizes are discrete random 
variables, so the arrival process is a (time-dependent) compound 
Poisson; see Feller (1967). Let B (u ) be the size of a bunch 
that arrives at time u; its distribution and generating function 
are, respectively, 

P{B (u) = k} = b k (u) 

oo 

E [ z B ^ ] = h ( z,u) = l 

k=l 



k = 1,2,... 



(3.1) 



z k b k (u) 



Again, let g(z,t,u) be the generating function of N(t,u) as 
in (2.7). Notice that g(z,t,u) represents the contribution 
to N(t,u) arising from arrivals in (u, u + du) plus the 
independent contribution from (u + du, t) ; see (2.8) . To 
derive the g.f . of the contribution from (u, du) , denoted by 
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AN, condition on the event of an arrival and bunch size B(u); 
by independence of the services. 



E [ z^ N | arrival in (u,du), B(u)] 

= [zF(t-u,u) + F(t-u,u)]^ U ^ 



Now removal of the condition on B(u) yields 

E [ z^ | arrival in (u, u + du) ] = h[(z-l) F(t-u,u) + l,u]. 



Of course 



„ r AN , , . , , . , 0 , 

E[z | no arrival in (u, u + du) ] = z = 1, 



so 



E [ z N ] = A (u) du h[( z-1) F(t-u,u) + 1 , u] + 1 - A (u) du + o(du) 



and thus 

g(z,t,u) = [A(u)du h[(z-l) F(t-u,u) +l,u] + 1 - A (u) du] g(z. 



+ o(du) , 

which leads to the differential equation 

- ^2- = A (u) (h[(z-l) F ( t-u , u) + l,u] - 1} g ( z , t , u) , 



(3.2) 

(3.3) 



(3.4) 

, u+du) 

(3.5) 

(3.6) 
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The solution of which is the generating function 



t 

g(z,t;u) = exp[- / A (v) {l-h[(z-l) F(t-v,v) + 1, v] }dv] (3.7) 

u 



Although inversion of (3.7) to produce a simple and familiar 
expression for the probabilities {p (t) , n = 0,1,2,...} seems 
beyond the realm of possibility, certain facts do emerge. 

For example, differentiation of g(z,t;0) yields moments, 



E[N(t) ] 




z=l 




X (v) E [ B ( v ) ] 



F(t-v,v)dv , 



(3.8) 



Var [N ( t ) ] = E [ N ( t ) ] + / X (v) E [B ( v) • (B (v ) -1 ) ] (F ( t-v , v) ) 2 dv, 

0 

and so forth. Note that Var [N ( t) ] /E [T] >1 if bunch sizes are 

sometimes greater than unity, as is to be expected. Furthermore, 
setting z = 0 yields 



t 

P{N (t) = 0} = p n (t,0) = exp{- / A (v ) h[F(t-v,v),v ]dv}(3.9) 

U 0 

The form of the general coefficient of z ^ may be deduced from 
consideration of (3.2) and (3.3). The coefficient of z- 1 in 
h[(z-l) F(t-u,u) + l,u] is the probability of exactly j 
"successes "--meaning survivals to time t from u in B(u) 
Bernoulli trials, hence 
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(3.10 



h[(z-l) F(t-u,u) + 1, u) ] = l z^c . (u) 

j=0 3 

00 k - k-s 

c . ( u) = l b,(u) (.) (F(t-U,u))- 1 (F ( t-u , u ) ) J , j = 0,1,2,. 

3 k=0 K 11 



Thus an alternative expression for g is 



g(z,t,-u) = exp[- / A(v) (1 - £ z^c.(v)}dv 

j=0 3 



U 



Since (c^ (v) , j = 0,1,2,...} is a discrete probability distri- 
bution for every v, as is evident from (3.10), N(t,u) clearly 
has a compound Poisson distribution. 

Of interest is the following 



EXAMPLE . Suppose A(v) = A , F(x,v) = e ^ x , b^(u) = (l-a)a^ ^ , 
k = 1,2,... . This is the case of stationary "stuttering Poisson" 

arrivals. Substitution into (3.7) yields for u = 0, 



g(z,t;0) = exp 



= exp 



(1-ct) [ (z-l)e' 9 (t v) +l] 
1 - a [ (z-l)e~° (t_v) +l] 



= exp 



[• C x l x 

r . f (z-l) e ~ 9 (t ~ v) dv 1 

|_ 0 (l-o- a[(z-l)e' S(t - v, J J 

1 r — a ~ a L ( ^t' |1 

( 1 - a - a ( z-l ) e ) J 



dv 



] 



= f 1 ~ - ct (z-l)e t 1 A / Ba 

|_ 1 - az J 
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as t tends to infinity the latter generating function approaches 



g ( z , °° ; 0 ) = 1 im g ( z , t ; 0 ) 

t -> 00 




A/6a 



the generating function of the negative binomial distribution. Hence 
the long-run distribution of server occupancy is, in this particular 
case, a familiar form that may readily be used in various applications 
in place of the ordinary Poisson that results from (2.6) or (2.10) 
under similar circumstances. 



4 . Bunch Division into Two Classes 

Suppose that bunches arrive in a Poisson manner, but that 
each bunch is independently fragmented into subbunches of type 1 
and type 2 customers with probability p^ and p^ respectively 
(although these probabilities may be time -dependent also, we do not 
oother with this) . Items of type i are served in accordance with 
distribution F.(x,u). This setup may model demands on certain 
logistics systems, e.g. by landing aircraft with different failure 
categories. Although the number of types of arrivals is limited 
to two, there is no difficulty in extending it to more types if 
necessary . 

Following the pattern leading to (3.5) one may write 
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] = g(z 1 ,z 2 ,t;u) 



N 1 (t,u) N 2 (t,u) 

1 Z 2 

= {X(u)du h[p 1 (z 1 -l) F^tt-U/U) +p 2 (z 2 ~l) F 2 (t-u,u) +1 ,u] +1-X (u) du} 

x g ( , z 2 , t ; u + du) (4.1) 

lb (t,u) is the contribution to system occupancy of the 
arrivals between u and t. It then follows that 

-X(u) {1 - h[p^(z^-l) F 1 (t-u,u) + p 2 (z 2 ~l) F 2 (t-u,u) + l,u] 

X g (Zj^ , Z 2 , t ; u) } , (4.2) 

from which the joint generating function 

g ( z ^ / z 2 ftj 0) 

t _ _ 

=exp[- / X (v) {1-h [p 1 (z 1 ~l) Fj^ (t-v, v) + p 2 (z 2 ~l) F 2 (t-v, v) + 1, v] }dv 

(4.3) 



where 
type i 

dg = 
du 



appears. Moments are obtained oy differentiation; see (3.8); 
to obtain means and variances simply replace X(u) by p^X(u). 

The covariance results from partial differentiations of the 
exponent at z^ = 1, z 2 = 1: 

t 

cov [Nj^ (t) ,N 2 (1;) ] = / x(v) E [B (v) (B (v) -1) ] P-]^ F^(t-v,v) F 9 (t-v,v)dv 

, . . . . (4.4) 

which is always positive. 
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EXAMPLE. Let A (u) = A, h(z,u) = (1-a) z/ (1-az) . If 
-0.x 

F^(x,u) = e then the following results: 



g(z 1 ,z 2 ,t;0) = exp 





-0 v -e V 


t ' 
-A / < 


p(z 1 ~l)e + p 2 (z 2 -l) e 


> 

< 

CD 

1 

> 

r 

CD 

1 


0 I 


[ 1 - a + a[p 1 (z 1 -l)e +p 2 (z 2 -l)e ] 



dv 



( 4 . 5 ) 



which, regrettably, cannot be integrated in closed form unless 
0^ = ©2 (=1 for convenience) in that case we find that as 

t ->■ CO 

] A/ a 

( 

___ 

This is the generating function of a bivariate distribution with 
negative binomial marginal distribution. 
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5 . Related Problems, or "Sons and Daughters of M/G/°°" 

Models for shot-noise, see Rice (1954) , and for dams and 
rainfall and runoff, see Gaver and Miller (1962>, share the general 
structure of the previous infinite server models. Time-dependent 
versions of these will be formulated and briefly discussed using 
the backward equation approach. 

Let A(u) be the rate of arrival of a certain event at time u, 
and let e(u,t) denote the random effect at time t of an event 

r 

at time u, 0 u t . In general e(u,t) will be real -valued 
random variable; in the shot noise application it represents the 
response at time t of an electrical circuit to an impuse at 
time u, and in the case of a dam or storage system it may be the 
amount of water in the reservoir resulting from a rainstorm at 
time u. Let the Laplace transform of e(u,t) be 



4> (S , t ;u) = E [e 



-S.§.(u, t) 



(5.1) 



Now write down a backward differential equation for 



ip (s , t;u) = E [e 



-sX ( t , u) 



] , 



(5.2) 



X(t,u) being the combined effect at t of all of the (sub, 
or component) effects occurring after u and before t: 



A(t ,u) 



t 

e(u.,t) = / e(s,t) dN(s) 



X ( t , ) = l 



(5.3) 



i=0 



u 
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where A(t) is the number of (Poisson) events in (u,t), and 
is the instant at which the ith such event occurs. It is seen that 



44 = [ A (u) - A (u) (f> ( s , t ; u) ] ip ( s , t ; u ) , (5.4) 

o U 

exactly as was true for (3.6) , and thus 

t 

ip(s,t;u) = exp { / A (u) [<p ( £ , t; v) -1 ] dv) (5.5) 

u 

EXAMPLE. Let A(u) = A and e,(u,t) = Se -6(t_u) , 0 > 0, 

S having the exponential distribution with density 
ye ,1X . It follows that 



<p ( s , t ; u ) 



y + se 



-0 (t-u) 



and thence that 



ip ( s , t ; 0 ) = 



( t 

exp / A 

( 0 



se 



-0 (t-u) 



y + se 



y + se 

-0 1 \ A/0 



y + s 



-e (t-u) du | 



A/9 



y + s 



(5.6) 



(5.7) 



as t so in the long run the total effect has gamma 

distribution. See Gaver and Miller (1962) for the same result 
derived differently. 
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5.1. 



Zipf's Law and Pareto Tails 

Let A(u) = e au and let e(u,t) = S(u) e 0 ^ , 6 > 0, 

S(u) having the distribution F( ) ; successive S(u) values are 
independent. This setup models a collection of organisms that 
are born at random times and grow independently and exponentially 
thereafter. We are interested in the fraction of all those born 
in (0,t) that exceed size x at time t; we shall see that the 
fraction exhibits the "Pareto tail" associated with Zipf's law; 
see Mandelbrodt (1978). 

Let 3!(x,u,t) denote the indicator function 

! 1 if e(u,t) >x 

(5.3) 

0 if e (u , t ) <_ x . 

In the present model define 

p (u,t) 5 E [ I (x , u , t ) ] = F (xe S ^ t ) (5.9) 

although what follows next does not require the latter explicit form. 
Now introduce the bivariate generating function 



g (z , z , t;u) 

d A 



E[z A(t ' u) 




(t ,u) 

] 



(5.10) 



where A(t,u) is the number of arrivals (births) in (u,t), and 
N (t,u) is the number of those organisms born in (u,t) that 
exceed x in size at time t; we put A(t,0) = A(t), and 
N (t) = N (t,0) . Then by the backward argument analogous to that 

X A 

producing (2.8) , 

g ( 2 a ' 2 x ' t •" u ) 

= z A(u)du p (u,t) + z A(u)du (1-p (u,t)) +1-A(u)du] (5.11) 

d ^ x a x 

X g (z , z , t;u + du) + o(du), 

a. A 
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which leads to a differential equation with solution 

t 

g(z a ,z x ,t;0) = exp {/ A (u) du [ z g ( z x ~l ) p x (u , t ) + (z a -l)]} (5.12) 

This shows that A(t) and N (t) have a bivariate Poisson distri- 

X 

bution; from (5.10) one finds 

t 

m(t) = E[A(t)] = Var [A (t ) ] = / A(u)du , 

0 

t 

m (t) = E[N (t) ] = Var [N (t ) ] = / A (u) p (u,t)du (5.13) 

XX X q X 

t 

Cov [A (t ) ,N (t) ] = / A (u) p (u,t)du = m (t) . 

x 0 x 



Under many interesting circumstances, a notable instance being 
the specific model beginning this example, both E[A(t)] and 
E [N (t) ] -* °° as t -*■ °°. 

X 



E [A (t ) ] = / e aU du = - (e at -l) ~ - e at , 

g a a 



E[N (t) ] = / e au du F(xe 
x 0 



-0 (t-u) 



= e at i f 

ex a/e -et 

xe 



a/0-1 = , > , 
z F (z) dz 



x /n , 

e r a/ 0—1 — , , 

~ tt- J z F(z)dz 



0x a/6 0 



provided the integral exists. 



(5.14) 
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Now define 



f ^ X,t) E [ A (t ) ] 



E[N (t) ] 

X 




F ( z ) dz , (t -> °°) (5 . 15 



the long-time average fraction of those organisms born before t 
and that exceed x in size at t; clearly for large x this 



the integral exists for large x, as will be assumed. 

It will now be shown that with high probability the above law 
should actually hold for observed data in the following sense. 

Form the ratio of observable random variables N (t)/A(t); this 

X 

ratio should approximate to f(x,°°) as t °°. To show that 
this is so, consider, for e > 0, 



For the specific model of this example both E[A(t)] and 
E [N (t) ] -*■ 00 as t -*• °°, consequently it can be shown (e.g. 
by the continuity theorem for characteristic functions) that 
(A(t),N (t) ) are approximately bivariate normal for large t 
with parameters given by (5.13). Therefore 





<_f(x,t) + e j =P{N x (t)-A(t)[f(x,t)+e] ^ 0 } . (5.1 



) 
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E[N (t) - A (t) [f (x, t) + e]] = - E[A(t)]e 

Var [N (t) - A ( t ) [ f (x , t ) + e]] 

X 

= Var [N (t) ] + [ f ( x , t ) + e] 2 Var[A(t)] 

X 

- 2[f(x,t) + e] covtN^Ct), A ( t ) ] 

= m (t) + [f(x,t) + e] 2 m(t) - 2[f(x,t) + e] m (t) 

X X 

= m(t) (f(x,t) + [f(x,t) + e] 2 - 2[f(x,t) + e] f(x,t)} 

= m(t) {f(x,t)(l - f(x,t)) + e 2 } (5.17) 



w use the normal approximation to assess the probability (5.16): 



P 



N (t) 

— < f (x,t) 

A ( t ) 



+ e = P{N (t) - A (t) [f (x,t) + e] <_ 0} 

) x 



^ 1 

/ 2 ¥ 



m(t)e/[m(t){f(x,t) (l-f(x,t)+e 2 }] 1//2 „ 2 /o 

/ e" z /2 dx 



(5.18 ) 



t -> 0 the fraction f(x,t) approaches the (finite) right side 
(5.15) and then, since m(t) ->• , the integral approaches unity; 

similar argument shows that N (t)/A(t) > f(x,t) - e with 

X 

obability approaching unity. It follows that the ratio 
(t)/A(t) is a consistent estimator of f(x,t) at the point x 
t cc. The statement (5.18) can also be used to supply 
proximate confidence limits for f(x,t). 
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An alternative formulation leads to similar asymptotic results. 
Suppose that the arrival rate is now taken to be kA (u) , k being 
a parameter that will later approach infinity; see Barbour (1974) 
for an analogous model and analysis. The interpretation is that 
when k becomes large organism births occur thick and fast — even 
more so, of course, for later times than earlier when (now) 

A (u) = ke as in our example. Now all analysis goes through as 
before, and the average fraction function is 



f k (x,t) = 



E[N (t) ] 



x 



E[A(t) ] 



1 

~a7T 



at 



x 

/ 



e<e at -n xe -0t 



z a/e F(z)dz , (5 



independent of k, while 



m ( t ; k ) = E [A ( t ) ] = k / A(u)du = - (e at -l) = km(t) 

0 a 



(5 



As k -*■ °° a central limit theorem argument once again applies 
(here for every finite t) to show that 



N (t) 

TTtT i £ k (x ' t) + e 



1 

/2tt 



km(t)e/[km(t) {f R (x,t) (l-f R (x,t) +e 2 ) ] 1/2 _ 2. 

/ e 7 dz 



(5 



and the latter probability clearly approaches unity as k + «, this 
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time for every t. The asymptotic normality also allows approximate 
confidence limits to be placed on f^(x,t) . 

We emphasize that the above analysis applies just to any single 
x-value. Analogous results should be derivable for any finite 
sequence of x-values, and thence extended by continuity to all 
real values, obtaining results similar to the Glivenko-Cantell i 
theorem for ordinary distribution functions. This, and other, 
generalizations are under development and will be reported in 
a subsequent paper. 
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DR. H. KOBAYASFI 

IB M 

YCFKTCWN HEIGHTS 
NEW YORK 


1 




10598 


CR. JOHN LEFOCZKY 
STATISTICS DEPART RENT 
CARNEG IE-MELLCN UNIVERSITY 
FITTS BURGH 
PENNSYLVANIA 


1 




15213 


L IERAR Y 
CCCE 55 

NAVAL POSTGRADUATE SCHOOL 
MONTEREY 


1 


CALIFORNIA 


93940 


DP. A. LEMOINE 
1020 GUINCA ST. 
PALO ALTO, 
CALIFGRN IA 


1 



943 01 



CR. J. MACCUEEN 
UNI V. OF CALI F. 
LOS ANGELES 
CALIFORNIA 


1 




90024 


FFCF. K. T. MARSFALL 
DEFT. OF OF 

NAVAL FOSTC-RACUATE SCHOOL 

MONTEREY 

CALI FOPNIA 


1 




53 940 


DR. M . MA2UMCAR 
MATH. DEPT. 
EST1NGF0LSE RES. LABS 
CHURCHILL BCFC 
FI TTSBLRGH 


1 


PENNSYLVANIA 


15235 


DP. LEON F. MCGINNIS 

SCHOOL OF INC. ANC SYS. ENG. 

GEORGIA INST. OF TECH. 

ATLANTA 

GEORGIA 


1 




30332 
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CR. 0. R. MCNEIL 
DEFT. CF STATISTICS 
PRINCETON UNIV. 

PRINCETON 
NEW JERSEY 

G8540 



PRCF. P. R. MILCH 
CEFT . CF OP 

NAVAL POSTGPACUATE SCFCCL 

MONTEREY 

CAL I FORN I A 

93940 



CR. F. MCSTELLER 
STAT. CEFT. 

FARVARC LNIV. 

CAMBRICC-E 

MASSACHUSETTS 

02139 



FPCF. R. R. READ 
DEFT . CF CP 

NAVAL FCSTGRACUATE SCHOOL 

MONTEREY 

CALI FORM I A 

93940 



DR . M . REI SER 
I B N 

THOMAS J. WAT S C N FES. CTR. 

YCRKTOWN HEIGHTS 
NEW YCFK 

10598 



DEAN CF RESEARCH 
CODE 0 13 

NAVAL FCSTGFACL ATE SCHOOL 
MONTEREY 

CALIFORNIA 93940 



FRCF . F. F. RICHARDS 
DEFT. CF OP 

NAVAL FCSTC-FACLATE SCHOOL 

MCNTERY 

CALIFORNIA 

93940 



DP. J. RICPCAN 
CEPT . OF MATHEMATICS 
ROCKEFELLER UNIV. 

NEW YORK 
NEW YORK 

100 21 



DP. LINUS SCHPAGE 
LNIV. CF CHICAGO 
GRAD. SCHOOL OF BLS. 
5836 GPEENWCCC AVE. 
CHICAGO, ILLINOIS 



60637 



No. of Copies 
1 



1 



1 



1 



1 



1 



1 



1 



1 
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No. of Copies 



DR. PAUL SCH l* 5 I T ZE R 

THCMAS J. WATSCN R cS EA F CH CTR 9 

PCST CFFICE BOX 218 

YORK TO NN HEIGHTS 

hE Vy YCFK 

10598 



CR. R ! CF £ F C S C R E N SCN 
CODE 303 NPRDC 
271 CATALINA BLVD. 

SAN DIEGC 
CALIFORNIA 

92152 



FFCF. M. G. SOVEREIGN 
DE FT . CF OF 

NAVAL FOSTER A CL AT E SCFCCL 

MONTEREY 

CALI FCRN IA 

53 94 0 



CR. V. SFINIVASAN 
GRADUATE SCHCCL CF BUSINESS 
STANFORD UN IV ERS IT Y 
STANFORD 
CALI FCRN I A 

94305 



DR . R . M. STARK 

STATISTICS ANC CCFPUTEF SCI. 

INIV. CF DELAWARE 

NE NARK 

CELAWARE 

19711 



FFOF. RICFARC VANSLYKE 
RES. ANALYSIS CCRP. 

6EECFWC0C 

CLC TAFFFN FOAC 

GLEN CCVEi NEW YORK 

11 542 



PRCF • JOHN N. TUKEY 
FINE HALL 
PRINCETON UNIV. 

PRINCETON 
NEW JERSEY 

08540 



CR . THCMAS C. VARLEY 
CFFICE OF NAVAL RESEARCH 
CODE 434 
ARLINGTON 
V A 

22217 



PRCF. G. H ATS ON 
FINE HALL 
PRINCETON UN IV . 

PRINCE TON 
NEN JERSEY 

• C8540 



1 

1 

1 

1 

1 

1 

1 

1 

1 
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